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ABSTRACT 

We present the results of a strong- lensing analysis of a complete sample of 12 very 
luminous X-ray clusters at z > 0.5 using HST/ACS images. Our modelling technique 
has uncovered some of the largest known critical curves outlined by many accurately- 
predicted sets of multiple images. The distribution of Einstein radii has a median value 
of 28^' (for a source redshift of Zs ~ 2), twice as large as other lower-z samples, and 
extends to 55" for MACS J0717. 5+3745, with an impressive enclosed Einstein mass 
of 7.4 X 10^^ Mq. We find that 9 clusters cover a very large area (> 2.5n') of high 
magnification (/i > xlO) for a source redshift of Zs ~ 8, providing primary targets 
for accessing the first stars and galaxies. We compare our results with theoretical 
predictions of the standard ACDM model which we show systematically fall short 
of our measured Einstein radii by a factor of 1.4, after accounting for the effect of 
lensing projection. Nevertheless, a revised analysis once arc redshifts become available, 
and similar analyses of larger samples, are needed in order to establish more precisely 
the level of discrepancy with ACDM predictions. 

Key words: dark matter, galaxies: clusters: individuals: MACS z > 0.5 sample, 
galaxies: clusters: general, gravitational lensing 



1 INTRODUCTION 

Massive galaxy clusters provide several means for indepen- 
dently examining any viable cosmological model. Cluster 
samples used for this purpose are usually complete in terms 
of X-ray and redshift measurements such as the RDCS 
(Rosati et al. 1998, Borgani et al. 1999) and the MACS 
survey (Ebeling, Edge, & Henry 2001, Ebeling et al. 2007) 
that we examine here with lensing. Ideally, controlled sam- 
ples of clusters with reliable lensing masses would allow the 
most direct comparison with theory but these are presently 
still in their infancy with only a handful of clusters discov- 
ered in wide-field weak lensing searches (Taylor et al. 2004, 
Wittman et al. 2006, Massey et al. 2007, Miyazaki et al. 
2007), and may be often subject to projection biases (e.g., 
Hennawi & Spergel 2005) . Weak lensing surveys should have 
the advantage of being approximately volume-limited over 
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the redshift range O.li^il.O, where the lensing distance ratio 
has a broad maximum. 

Ongoing measurements of the Sunyaev-Zel'dovich (SZ) 
effect (which has already been detected in many tens of clus- 
ters, including some at relatively high redshifts z ^ 0.9, 
e.g., LaRoque et al. 2006) in a large number of clusters over 
a wide redshift range should improve upon X-ray selection, 
due to the redshift-independence of the effect (though clearly 
high redshift clusters are relatively less well resolved). Also, 
since the insensitivity of the effect to internal gas density 
variations will not be biased as with X-ray flux selection 
towards systems with luminous shocked gas, measurements 
should better correlate with cluster masses. 

Strong gravitational lensing (SL) is nearly always seen 
in suflficiently detailed observations of massive clusters. Ein- 
stein radii derived from such observations provide a reliable 
projected mass within the observed Einstein radius, which 
depends only on fundamental constants, G and c, and knowl- 
edge of the lens and source distances. Constraining the inner 
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mass profile from SL is mucii iiarder, requiring many sets of 
multiple images covering a wide range of source redshift, to 
overcome inherent lensing degeneracy between the gradient 
of the mass profile and the scaling of the bend angle with 
source distance, so that SL based mass profiles have been 
usefully constrained for only several clusters producing reli- 
able mass maps (e.g., Abell 370, Kneib et al. 1993, Richard 
et al. 2009a; Abell 901, Deb et al. 2010; Abell 1689, Broad- 
hurst et al. 2005, Coe et al. 2010; Abell 1703, Limousin et 
al. 2008; C10024+1654, Liesenborgs et al. 2008, Zitrin et al. 
2009b; MS 2137.3-2353, Gavazzi et al. 2003, Merten et al. 
2009; RXJ1347, Bradac et al. 2008a, Halkola et al. 2008; 
SDSS J1004+4112, Sharon et al. 2005, Liesenborgs et al. 
2009; "The bullet cluster", Bradac et al. 2006). 

In quite a few of these clusters there is a significant in- 
consistency between the Einstein radius directly measured 
from SL analyses and the Einstein radius predicted from 
model profiles fitted to weak-lensing measurements of the 
same clusters, often subject to background selection and di- 
lution problems over a wide radial range (e.g., Medezinski et 
al. 2010), though more so towards the SL regime where the 
degree of contamination by cluster members is often poorly 
corrected. The importance of a SL analysis of a significant 
sample is clear, especially in the statistical aspect, where 
the modelling method is similar for all examined clusters 
and thus supplies a coherent view for objective comparison 
(e.g., Richard et al. 2009b). 

The development of SL modelling-methods has in- 
creased in response to the improvement in data quality since 
the discoveries of the first giant arcs (e.g., Lynds & Petrosian 
1986, Soucail et al. 1987, Kneib et al. 1996). Most methods 
can be classified as "parametric" if based on physical param- 
eterisation, and as "non-parametric" if they are "grid-based" 
(see also §4.4 in Coe et al. 2008, and references therein). Still, 
many of these methods include too many parameters to be 
well-constrained by the number of initially known multiply- 
lensed systems. Here we use ACS imaging to identify new 
multiply-lensed systems in order to constrain the mass dis- 
tributions and define the critical curves of the sample, moti- 
vated by the successful minimalistic approach of Broadhurst 
et al. (2005) to lens modelling. In following work (Zitrin et al. 
2009b) we presented an improved modelling method which 
involved only 6 free parameters enabling easier constraint, 
since the number of constraints has to be equal or larger to 
the number of parameters in order to get a reliable fit. Two 
of these parameters are primarily set to reasonable values 
so only 4 of these parameters have to be constrained ini- 
tially, which sets a very reliable starting-point using obvious 
systems. The mass distribution is therefore primarily well 
constrained, uncovering many multiple-images which can be 
then iteratively incorporated into the model, by using their 
redshift estimation and location in the image-plane. 

The Massive Cluster Survey (MACS) has been aimed to 
create a complete sample of the very X-ray luminous clusters 
in the Universe, successfully increasing the number of known 
such clusters to hundreds, particularly at z > 0.3 (Ebeling, 
Edge, & Henry 2001). From this a complete sample of 12 
high-z MACS clusters {z > 0.5) was defined by Ebeling 
et al. (2007) which have proved very interesting in several 
follow-up studies including deep X-ray, SZ and HST imag- 
ing. Detection of a large-scale filament has been reported in 
the case of MACS J0717.5-h3745 by Ebeling et al. (2004), 



for which many multiply-lensed images have been recently 
identified by Zitrin et al. (2009a), revealing this object to 
be the largest known lens, with an Einstein radius equiva- 
lent to 55'' (for a source at z ~ 2.5). In the case of MACS 
J1149. 5+2223 (Zitrin k Broadhurst 2009), a background 
spiral galaxy at z = 1.49 (Smith et al. 2009) has been shown 
to be multiply-lensed into several very large images, requir- 
ing a very shallow, unrelaxed central mass distribution. An- 
other large multiply-lensed sub-mm source at a redshift of 
z ^ 2.9 has been identified in MACS J0454. 1-0300 (also re- 
ferred to as MS 0451.6-0305; Takata et al. 2003, Borys et 
al. 2004, Berciano Alba et al. 2007, 2009), MACS J0025.4- 
1222 was found to be a "bullet cluster" -like (Bradac et al. 
2008b), and other MACS clusters have been recently used 
for an extensive arc statistics study (Horesh et al. 2010). 

The X-ray data available for this sample (see Ebel- 
ing et al. 2007) along with the high-resolution HST/ ACS 
imaging and SZ imaging (e.g., Laroque et al. 2003) make 
these MACS targets particularly useful for understanding 
the nature of the most massive clusters. Here we complete a 
full SL analysis of this 12 cluster sample with the available 
deep HST/ ACS imaging, principally to derive their Einstein 
radii and also to help constrain the central mass distribu- 
tions. The effective Einstein radii derived here are the corre- 
sponding radii of circles of equivalent areas to those enclosed 
within the critical curves, or similarly, the radii within which 
the averaged k, is equal to 1. In addition, this work can also 
supply detailed magnification maps to help motivate deeper 
searches for high-z galaxies behind these clusters. 

Another major motivation for pursuing accurate 
lensing-maps is the increased precision of model predictions 
for cluster- size massive halos in the standard ACDM model 
and close variance (see Umetsu & Broadhurst 2008, e.g.. 
Bullock et al. 2001, Spergel et al. 2003, 2007, Tegmark et 
al. 2004, Hennawi et al. 2007, Neto et al. 2007, Duffy et 
al. 2008, Keselman, Nusser k Peebles 2009, Meneghetti et 
al. 2010, Fedeli et al. 2010). In the standard hierarchical 
model, large massive bodies are due to collapse recently 
and thus are not expected to be found in large numbers 
at high redshifts. On the other hand, the larger volume 
available with increasing distance means that in practice 
we cannot expect to reside next to the most massive cluster 
(see also Zitrin et al. 2009a). The increasing number of ac- 
curately analysed clusters have revealed larger lenses than 
predicted by the standard ACDM model (Broadhurst k 
Barkana 2008, Puchwein k Hilbert 2009). This discrepancy 
is empirically supported by the surprisingly concentrated 
mass profiles measured for such clusters, when combining 
the inner strong lensing with the outer weak lensing sig- 
nal (Gavazzi et al. 2003, Broadhurst et al. 2005 k 2008, 
Limousin et al. 2008, Donnarumma et al. 2009, Oguri et 
al. 2009, Umetsu et al. 2009, Zitrin et al. 2009b,2010), or 
independently, when using the internal dynamics of cluster 
members (Lemze et al. 2009) and deep X-ray data (Lemze 
et al. 2008). Geometrically, lensing is considered to be op- 
timised at intermediate redshifts, where for a given mass 
the critical density for lensing is minimal, but this is partly 
offset by the late hierarchical growth of high-mass systems. 
This trade-off results in estimates of the amplitude of strong 
lensing to favour the redshift range z = 0.2 — 0.4. Our recent 
analysis of two X-ray selected MACS clusters reveals very 
large lenses (e.g., MACS J1149.5-h2223, Zitrin k Broadhurst 
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2009; MACS 0717.5+3745, Zitrin et al. 2009a) at high red- 
shift, increasing the tension with ACDM predictions, since 
large Einstein radius clusters at high- 2; require earlier for- 
mation than implied by the standard ACDM model (Sadeh 
& Rephaeli 2008). The largest lensing clusters have proven 
to be excellent targets for accessing the faint early Universe 
due to their large magnification consistently providing the 
highest redshift galaxies (Ebbels et al. 1996, Franx et al. 
1997, Frye & Broadhurst 1998, Bouwens et al. 2004, Kneib 
et al. 2004, Bradley et al. 2008, Zheng et al. 2009). 

The statistics of large Einstein radii provide an impor- 
tant opportunity to test the standard ACDM paradigm, as it 
probes both the high-mass end of the cluster mass function 
and central mass distributions of massive clusters (Oguri & 
Blandford 2009). Increasing numbers of theoretical and ob- 
servational works have been done recently, deriving Einstein 
radius distributions of various samples. For example, Cypri- 
ano et al. (2003) have derived the mass distributions of 24 
X-ray selected Abell clusters via weak-lensing. By fitting to 
SIS profiles they obtain a mean Einstein radius of 17" . 
More recently, Hoekstra (2007) finds by weak-lensing anal- 
ysis and SIS fitting a mean Einstein radius of ^ 14" for a 
sample of 20 X-ray luminous clusters. Okabe et al. (2009) 
find the mass distribution and Einstein radii for a sample 
of ^ 30 LoCuSS clusters by fitting the weak-lensing data to 
CIS profiles, obtaining a mean value of ^ 22^^ (see Okabe 
et al. 2009 and references therein). Note that these values 
quoted here are calculated from the corresponding tables in 
these papers, and are not necessarily scaled or normalised 
to certain lens and source redshift s. 

Other recent work more relevant to our study sum- 
marises SL analysis for a sample of 20, mostly relaxed undis- 
turbed clusters (Richard et al. 2009b). They find that the 
Einstein radii are distributed log-normally with a peak at 
6e — 14.45" and a corresponding 1.95 x lO^^M© enclosed 
mass (within R<250 kpc), and show that the predicted dis- 
tribution of Einstein radii from ACDM cosmology falls short 
of the observed Einstein radii by a factor of 2. It has been 
claimed that much larger Einstein radii can be contemplated 
only with mass distributions which are highly prolate and 
aligned along the line of sight (Corless & King 2007, Oguri 
& Blandford 2009, see also Hennawi et al. 2007, Meneghetti 
et al. 2007, Sereno, Jetzer & Lubini 2010). Analysing the 
12 X-ray selected, high-^ MACS clusters is important since 
they are predicted to be massive and should form efficient 
lenses, whose properties can be then compared to such stud- 
ies, playing an important role in probing the ACDM sce- 
nario. 

The paper is organised as follows: in §2 we describe the 
sample and observations; the lensing analysis is described 
in §3 detailing each cluster separately. Our overall results 
are presented and discussed in §4, followed by a Conclusion 
(§5). Throughout the paper we adopt the standard cosmol- 
ogy {Qmo = 0.3, Qao = 0.7,/i = 0.7). 



2 THE SAMPLE AND OBSERVATIONS 

The z > 0.5 MACS clusters have been imaged with the 
Wide Field Channel (WFC) of the ACS installed on HST, 
mainly in the framework of proposal ID 9722 (PI: Ebeling). 
This mainly comprises two-filter observations in the F555W 



and F814W bands, taken during 2003 and 2004 with typi- 
cal exposure times of ^ 4500s in each filter. We make use 
of some other HST data available in the Hubble Legacy 
Archive (HLA), which we list in Table [l] In addition, avail- 
able SExtractor (Bertin & Arnouts 1996) photometry cata- 
logues were also downloaded from the HLA, and were used 
to construct colour-magnitude diagrams for identifying the 
red cluster sequence galaxies belonging to each cluster, as a 
starting point for each lensing model, as detailed in ^ 

3 STRONG LENSING MODELLING AND 
ANALYSIS 

We apply our well tested approach to lens modelling, which 
we have applied successfully before to various clusters un- 
covering unprecedentedly large numbers of multiply-lensed 
images in A1689, C10024, MACS J1149. 5+2223 and MACS 
J071 7. 5+3745 (respectively, Broadhurst et al. 2005, Zitrin 
et al. 2009b, Zitrin & Broadhurst 2009, Zitrin et al. 2009a). 
The full details of this approach can be found in these ear- 
lier papers. Briefly, the basic assumption adopted is that 
mass approximately traces light, so that the photometry of 
the red cluster member galaxies is the starting point for 
our model. In a recent paper (Zitrin et al. 2010) we anal- 
yse Abell 1703 and show that this assumption is well based, 
by comparison to an assumption-free non-parametric tech- 
nique (Liesenborgs et al. 2006) which yields a very similar 
overall mass distribution. In addition, we show that similar 
to other clusters we have analysed to date, our modelling 
method has the inherent flexibility to find and reproduce 
many multiple-images even if initially constrained by only a 
few obvious systems. 

As mentioned, the starting point of the model are clus- 
ter member galaxies, which are identified as lying close to 
the cluster sequence by the photometry provided in the Hub- 
ble Legacy Archive (in this process sometimes massive fore- 
ground galaxies are included as well, scaled down by the 
relative distance ratio, since they can locally affect nearby 
images) . We then approximate the large scale distribution of 
matter by assigning a power-law mass profile to each galaxy, 
the sum of which is then smoothed. The degree of smooth- 
ing (S) and the index of the power-law (q) are the most 
important free parameters determining the mass profile. A 
worthwhile improvement in fitting the location of the lensed 
images is generally found by expanding to first order the 
gravitational potential of the smooth component, equivalent 
to a coherent shear describing the overall matter ellipticity, 
where the direction of the shear and its amplitude are free, 
allowing for some flexibility in the relation between the dis- 
tribution of DM and the distribution of galaxies which can- 
not be expected to trace each other in detail. This freedom 
also allows the effective centre to slightly shift, as was the 
case in our analysis of C10024 (Zitrin et al. 2009b). The to- 
tal deflection field ariO), consists of the galaxy component, 
^gai{0), scaled by a factor Kgai, the cluster DM component 
o^dm{0), scaled by (l-Kgai), and the external shear compo- 
nent aex{0): 

aT0) = Kgaiagai{0) + (1 - Kgai)aDM{0) + (1) 

where the deflection field at position Om due to the external 
shear, aex(Om) = {oLex,x,o^ex,y), is given by: 
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Table 1. Properties of the MACS z > 0.5 sample. The following data (columns 1-7) are based on Ebeling et al. (2007): Column 1: cluster 
name in the MACS survey; Columns 2 & 3: The RA and Declination of the X-ray centroids (as determined from Chandra ACIS-I data); 
Column 4-' redshifts; Column 5: velocity dispersions, in km s~^; Column 6: Chandra X-ray luminosities in 10^^ erg quoted for the 
0.1-2.4 keV band. These luminosities are quoted within r2oo and exclude X-ray point sources; Column 7: X-ray temperatures, measured 
from Chandra data within riooo, but excluding a central region of 70 kpc radius around the listed X-ray centroid; Column 8: Ebeling 
et al. (2007) morphology code - assessed visually based on the appearance of the X-ray contours and the goodness of the optical/X-ray 
alignment. The assigned codes (from apparently relaxed to extremely disturbed) are 1 (pronounced cool core, perfect alignment of X- 
ray peak and single cD galaxy), 2 (good optical/X-ray alignment, concentric contours), 3 (nonconcentric contours, obvious small-scale 
substructure), and 4 (poor optical/X-ray alignment, multiple peaks, no cD galaxy), errors are estimated as less than 1; Column 9: ACS 
bands used here. Note, the clusters MACS J0018.5+1626 and MACS J0454.1-0300 are better known from earlier work as C10016+1609 
and MS 0451.6-0305, respectively. For more information see Ebeling et al. (2007). 



MACS a (J2000.0) 5 (J2000.0) z a L^^chandra KT{kev) M.C.E. ACS bands 



J0018.5+1626 


00 


18 


33. 


,835 


+16 26 16.64 


0, 


,545 




19.6 + 0.3 


9.4 + 1.3 


3 


F606W,F775W,F850LP 


J0025.4-1222 


00 


25 


29. 


,381 


-12 22 37.06 


0, 


.584 


740+90 


8.8 + 0.2 


7.1 + 0.7 


3 


F450W,F555W,F814W 


J0257.1-2325 


02 


57 


09. 


.151 


-23 26 05.83 


0, 


.505 




13.7 + 0.3 


10.5 + 1.0 


2 


F555W,F814W 


J0454. 1-0300 


04 


54 


11, 


.125 


-03 00 53.77 


0, 


.538 


1250ilfo 


16.8 + 0.6 


7.5 + 1.0 


2 


F555W,F775W,F814W 


J0647.7+7015 


06 


47 


50, 


.469 


+70 14 54.95 


0, 


.591 


900+120 


15.9 + 0.4 


11.5 + 1.0 


2 


F555W,F814W 


J0717.5+3745 


07 


17 


30, 


.927 


+37 45 29.74 


0, 


.546 




24.6 + 0.3 


11.6 + 0.5 


4 


F555W,F606W,F814W 


J0744.8+3927 


07 


44 


52, 


.470 


+39 27 27.34 


0, 


.698 


instill 


22.9 + 0.6 


8.1 + 0.6 


2 


F555W,F814W 


J0911.2-hl746 


09 


11 


11, 


.277 


+17 46 31.94 


0, 


.505 




7.8 + 0.3 


8.8 + 0.7 


4 


F555W,F814W 


J1149.5-h2223 


11 


49 


35, 


.093 


+22 24 10.94 


0, 


.544 


1840tlfo 


17.6 + 0.4 


9.1 + 0.7 


4 


F555W,F814W 


J1423.8+2404 


14 


23 


47, 


.663 


+24 04 40.14 


0, 


.543 


1300+^7° 


16.5 + 0.7 


7.0 + 0.8 


1 


F555W,F814W 


J2129.4-0741 


21 


29 


26, 


.214 


-07 41 26.22 


0, 


.589 




15.7 + 0.4 


8.1 + 0.7 


3 


F555W,F814W 


J2214.9-1359 


22 


14 


57, 


.415 


-14 00 10.78 


0, 


.503 


1300tfoo 


14.1 + 0.3 


8.8 + 0.7 


2 


F555W,F814W 



OLex,x{0m) = ItI cos(2(/)^) Aa^^ri + ItI sin(2(/)^) Ay^ , (2) 

OLex,y{0m) = ItI Sm{2(f>^)AXm " ItI COs(20^) A^^ri , (3) 

and (Axm, A^rn) is the displacement vector of the position 
Om with respect to a fiducial reference position, which we 
take as the lower-left pixel position (1, 1), and 0^ is the po- 
sition angle of the spin- 2 external gravitational shear mea- 
sured anti-clockwise from the x-axis. The normalisation of 
the model and the relative scaling of the smooth DM compo- 
nent versus the galaxy contribution brings the total number 
of free parameters in the model to 6. 

Note that since our goal is to find the Einstein radius 
distribution and the masses within the critical curves, and 
since many clusters lack multiple-images redshift informa- 
tion, we do not attempt to accurately constrain here the 
mass profiles of the sample clusters and thus the parame- 
ters q and S are relatively irrelevant. Any reasonable values 
for these parameters yield similar critical curves, Einstein 
radii and thus also the corresponding masses enclosed within 
them. This effectively leaves us with 4 parameters per model 
which can be easily constrained by the multiple-images in- 
corporated and found here. Still, we have explored also the 
q and S parameter space for each cluster to verify this. In 
addition, we showed that in cases where there are too few 
bands in order to obtain reliable photometric redshifts and 
where no spectroscopic redshifts are available, a very effi- 
cient assumption is that the outer, blue multiple-images are 
at a redshift of 2: ~ 2 — 2.5. This assumption is based on 
previous analysis of many clusters for which the vast major- 
ity of the outer blue images are at this redshift range due 
to the nesting effect by which the critical curve expands for 
higher redshift objects, yet the redshift of these blue im- 
ages is limited by the Lyman-alpha break to be lower than 



z ^ 3. For example, this assumption was proven to be very 
accurate in MACS Jl 149.5+2223, where the intermediate- 
distance blue images we assumed at z ^ 2 (system 3 in Zitrin 
& Broadhurst 2009) were later spectroscopically found to 
be at z = 1.89, and our corresponding estimation of the 
spiral galaxy at z ^ 1.5 (system 1 in Zitrin & Broadhurst 
2009), was later verified spectroscopically to be at z = 1.49 
(Smith et al. 2009). In addition, in our works on A1689 and 
C10024 we showed that the same minimum is obtained both 
when minimising according to the images location in the 
image-plane, and when minimising according to the photo- 
metric redshifts. This means that in practice also the slope 
which generally can be constrained only using redshifts in- 
formation, can be well approximated by minimising solely 
the image-plane RMS of the reproduced images, especially 
if a clear minimum is seen (see Zitrin et al. 2009b). Note, 
here we aim to constrain the critical curves and the mass 
enclosed within them which are both relatively independent 
and indifferent to the profile, enabling an accurate measure- 
ment also without redshift information, which in turn can 
be used to approximate the magnification of the lens, since 
this quantity is profile dependent. 

Firstly we use our preliminary models to lens various 
well detected candidate lensed galaxies back to the source 
plane using the derived deflection fleld, and then relens this 
source plane to predict the detailed appearance and location 
of additional counter images, which may then be identified 
in the data by morphology, internal structure and colour. 
We stress that multiple images found this way must be ac- 
curately reproduced by our model and are not simply eyeball 
"candidates" requiring redshift verification. Note, due to the 
volume of this work we do not attempt to exhaust each clus- 
ter field finding great numbers of multiple-images, and con- 
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centrate instead on finding clear examples of multiply-lensed 
systems to constrain the models and obtain the Einstein ra- 
dius. The fit is assessed by the RMS uncertainty in the image 
plane: 

RMS images = ^((^z " + (l/z " 2/0^) / J^images, (4) 

i 

where and 1/^ are the locations given by the model, and 
Xi and i/i are the real images location, and the sum is over 
all Nimages imagcs. The best-fit solution is obtained by the 
minimum RMS, and the uncertainties are determined by the 
location of predicted images in the image plane. 

Importantly, this image-plane minimisation does not 
suffer from the well known bias involved with source plane 
minimisation, where solutions are biased by minimal scatter 
towards shallow mass profiles with correspondingly higher 
magnification. The model is successively refined as addi- 
tional sets of multiple images are identified and then incor- 
porated to improve the fit, using also their redshifts mea- 
surements or estimates for better constraining the mass 
slope through the cosmological relation of the Dis/Ds 
growth (see Zitrin et al. 2009b). We detail this procedure 
for each cluster in the corresponding subsection below, each 
followed by corresponding figures of the multiple-images and 
the critical curves, and the mass distribution (Figures [l] to 
24|. 



3.1 MACS J0018.5+1626 

The galaxy cluster MACS J0018. 5+1626 [z = 0.546, also 
known as CI 0016+1609, or MS 0016) has been subject 
to extensive study, mainly due to its high X-ray luminos- 
ity, strong SZ effect, and high redshift (e.g.. Dressier & 
Gunn 1992, Dressier et al. 1999, Luppino et al. 1999, El- 
lis & Jones 2002, Clowe et al. 2000, 2003, Laroque et al. 
2003), setting an example of a rich distant cluster (Worrall 
& Birkinshaw, 2003, and references therein). It has been es- 
tablished that this cluster is the major part of a large-scale 
filamentary structure (Worrall & Birkinshaw, 2003, and ref- 
erences therein, Tanaka, 2007). This fact and the promi- 
nent X-ray emission imply a high cluster mass, which should 
be manifested in multiply-lensed systems spread through- 
out the field. However, due to lower quality optical data 
and the complexity of its structure, no multiply-lensed sys- 
tems were found in MACS J0018. 5+1626 up to date. Only 
one strongly-lensed arc was referred to before in this clus- 
ter (Lavery 1996), yet no spectroscopic redshift nor counter- 
images were introduced. We begin our analysis by construct- 
ing an initial model in the method described above (^, 
where as a starting point we insert averaged values to the 
various parameters, since no system is known a-priori to 
constrain the fit. We then notice the similar colours and 
symmetry of the images of systems 1, 2, and 3. After ver- 
ifying that the reproduction of these systems is physically 
likely in the context of our model, we use these systems to 
fully constrain the model and derive the critical curves (see 
Figure [l]), their corresponding Einstein radius {Be — 28 + 2^^ 
for a source at ^ 2; systems 2 & 3) and the mass en- 
closed within them (1.46 + 0.1 x lO^'^M©), see Figures[l]and 
|2]for more information. Accordingly, the estimated redshift 
of system 1 is ^ 1.5. 




Figure 1. Galaxy cluster MACS J0018. 5+1626 {z = 0.546, also 
known as CI 0016+1609, or MS 0016) imaged with Hubble/ACS 
F606W, F775W, and F850LP bands. The blue curve overlaid 
shows the tangential critical curve corresponding to the distance 
of system 1 at an estimated redshift of z ~ 1.5, and which passes 
through the close triplet of lensed arcs in this system. The larger 
critical curve overlaid in white corresponds to the larger source 
distance estimated as z ~ 2 — 2.5, passing through the close pairs 
of the candidate images of systems 2 and 3. This critical curve 
encloses a large lensed region, with an equivalent Einstein radius 
of ~ 150 kpc at the redshift of the cluster. Note, our model pre- 
dicts tiny images of pairs 2.1/2.2 and 3.1/3.2 at the other side of 
the cluster, which we are not able to securely detect due to poor 
depth and lack of colour range. 

There are various other mass estimates for this cluster, 
mainly from weak-lensing analyses. Small et al. (1997) quote 
values of - 2.9 x lO^^M© and - 1.87 x lO^^M© within 200 
kpc (depending on the method, see also Small et al. 1995), 
and Hoekstra (2007) quotes a value of 7.9x lO^^M© within 
500 kpc, all in general accordance with our measurement. 
The Einstein radius for this cluster was also mentioned be- 
fore. Hoekstra (2007) quotes ^ 9'' from a SIS model fit to 
the tangential distortion from 0.25 to 1.5 Mpc, while 
Williams, Navarro & Bartelmann (1999) quote 25^^ accord- 
ing to the arc discussed by Lavery (1996), in accordance 
with our estimation. The reference centre of our analysis is 
fixed near the centre of the F775W (program ID 10493) ACS 
frame at: RA = 00:18:33.41, Dec = +16:26:14.95 (J2000.0), 
where one arcsecond corresponds to 6.42 kpc at the redshift 
of this cluster. 



3.2 MACS J0025.4-1222 

The galaxy cluster MACS J0025.4-1222 {z = 0.584) and 
its merging properties were recently analysed by Bradac et 
al. (2008b), which found multiply-lensed images of 4 back- 
ground galaxies, and obtained spectroscopic redshifts for a 
pair of them and photometric redshifts for the other two sys- 
tems: system A+B (system 1 here) at Zspec = 2.38, system 



6 Zitrin et al 





Figure 2. 2D surface mass distribution (k), in units of the crit- 
ical density, of MACS J0018. 5+1626. Contours are shown in hn- 
ear units, derived from the mass model constrained using the 
multiply- lensed images seen in Figure^ Axes are in ACS pixels 
(0.05" /pixel), and a 20" scale bar is overplotted. 



C (system 2 here) at 



■'phot 



1+0.5 



and system D (system 



3 here) at Zphot = 2.8l5"8- systems 1 and 2 (see 

Figure [3} and one more uncovered arc (system 4 in Figure 
[s] Zs ^ 1-9), to derive the corresponding Einstein radius 



30 di 2^' for Zs ^ 2.38) and the mass enclosed within 
it (2.421°; ^3 X lO^^M©), in excellent agreement with the re- 
sult of Bradac et al (2008b; 2.5 x lO^^M©, within 300 kpc 
centred on the SE BOG). In addition we find that system 
A in Bradac et al (2008b; system 1 here) includes also an 
additional arc on the other side of the BOG, as seen in Fig- 
ure[3](arc 1.2), which our model reproduces very accurately 
(see Figure 30). The reference centre of our analysis is fixed 
near the centre of the ACS frame at: RA = 00:25:29.61, Dec 
= -12:22:52.89 (J2000.0), where one arcsecond corresponds 
to 6.64 kpc at the redshift of this cluster. 



3.3 MACS J0257.1-2325 

In the galaxy cluster MACS J0257. 1-2325 (z = 0.505) vari- 
ous arcs are seen, from which we use 13 multiply-lensed im- 
ages, belonging to ~ 7 background distant galaxies, in the 
redshift range z ^ 1 to z ^ 4 to fully constrain the model. 
The reference centre of our analysis is fixed near the centre 
of the ACS frame at: RA = 02:57:07.24, Dec = -23:26:02.88 
(J2000.0), where one arcsecond corresponds to 6.17 kpc at 
the redshift of this cluster. This cluster was not analysed 
before and we did not find any arc redshift information. We 
assume a redshift of ^ 2 — 2.5 for the faint distant blue 
images (system 4 in Figure [s]) according to which the spec- 
tacular images next to the cD galaxy (system 1; see also 
Figure [3T| are at Zs I, the blue images of systems 2 and 
3 are at Zs 1.5 — 2, and the red drop-out candidates are 
at Zs ^ 3.5 — 4 (system 5) helping to constrain the fit. We 
find that the critical curves (Figure [s]) have an equivalent 
Einstein radius of 39 =b 2" (for Zs ^ 2) and enclose a mass 
of 3.35t°;?g X lO^^M©. . 




Figure 3. Galaxy cluster MACS J0025. 4-1222 (z = 0.584) im- 
aged with Hubble/ ACS F555W and F814W bands. The overlaid 
critical curve (white) corresponds to system 1 (zs = 2.38), enclos- 
ing a critical area of an effective Einstein radius of ~ 200 kpc at 
the redshift of this cluster. 




Figure 4. 2D surface mass distribution (k), in units of the crit- 
ical density, of MACS J0025. 4-1222. Contours are shown in hn- 
ear units, derived from the mass model constrained using the 
multiply-lensed images seen in Figure [3] Axes are in ACS pixels 
(0.05" /pixel), and a 20^^ scale bar is overplotted. 



3.4 MACS J0454.1-0300 

The galaxy cluster MACS J0454. 1-0300 (z = 0.538, also 
known as MS 0451.6-0305) has been subject to extensive 
study due to its high X-ray luminosity, strong SZ effect, 
and high redshift (e.g., EUingson et al. 1998, Molnar et al. 
2002, Laroque et al. 2003, Geach et al. 2006), and is known 
to host various lensed sub-mm sources with radio counter- 
parts (Takata et al. 2003, Borys et al. 2004, Berciano Alba 
et al. 2007, 2009), which we incorporate in our analysis be- 
low. Previous strong-lensing models of this cluster were in- 
troduced by Borys et al. (2004) and Berciano Alba et al. 
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Figure 5. Galaxy cluster MACS J0257.1-2325 {z = 0.505) im- 
aged with Hubble/ACS F555W and F814W bands. The tangen- 
tial critical curves overlaid in blue correspond to a source redshift 
of Zs ~ 1, passing through the close pair of images 1.2 and 1.3. 
The larger critical curves overlaid in white correspond to a source 
redshift of ~ 2 — 2.5, passing through the close pair of images 
4.1 and 4.2. 




500 1000 1500 2000 2500 

Pixels 



Figure 6. 2D surface mass distribution (^c), in units of the crit- 
ical density, of MACS J0257. 1-2325. Contours are shown in lin- 
ear units, derived from the mass model constrained using the 
multiply- lensed images seen in Figure [s] Axes are in ACS pixels 
{Q.^h" /pixel) ^ and a 20^^ scale bar is overplotted. 



(2009) with similar symmetry as our model. We fully con- 
strain our model with the arcs and their redshift information 
listed in these papers (see also Figure |7| here) , from which 
we derive an effective Einstein radius of 19 ± 2" for a source 
redshift of Zs — 2.9 (system 1 in Figure [7| and a mass of 
0.82lo oi X lO^^M© enclosed within this critical curve. For a 
source redshift of ~ 2 we correspondingly get an effective 
Einstein radius of 13 ±2'', and a mass of 0.4ll° °^ x lO^^M© 
enclosed within the critical curve (see Figure It]). Berciano 
Alba et al. (2009) quote a mass of 1.73 x lO^M© within 
30 " of the cluster centre. Our model yields 1.8 x lO^^M© 




Figure 7. Galaxy cluster MACS J0454. 1-0300 {z = 0.538, also 
known as MS 0451.6-0305) imaged with Hubble/ACS F555W, 
F775W, and F814W bands. The tangential critical curves overlaid 
in white correspond to a source redshift of Zg = 2.9 (system 1, see 
also Borys et al. 2004), passing through the close pair of images 
1.1/1.2, and 4.1/4.2. The smaller critical curves overlaid in blue 
correspond to a source redshift oi Zs ~ 2. 



within 30 " of the cluster centre, in full agreement. In our 
modelling process we find two other systems of multiply- 
lensed images. The first (system 3) is a triplet of faint arcs 
next to the core, at an estimated redshift oi Zs ^ 1.5 — 2. 
The second system (number 4) consists of 4 bright images, 
all have prominent emission in the K'-band (see Takata et 
al. 2003), corresponding to a redshift of Zs — 2.9 similar to 
system 1. See Figure |7| for more details. The reference cen- 
tre of our analysis is fixed near the centre of the F555W 
(proposal ID 9722) ACS frame at: RA = 04:54:10.36, Dec = 
-03:01:03.02 (J2000.0), where one arcsecond corresponds to 
6.37 kpc at the redshift of this cluster. 



3.5 MACS J0647.7+7015 

In the galaxy cluster MACS J0647.7+7015 {z = 0.591) we 
note a remarkable blue background system lensed 6 times, 
where each lensed image is a quartet of arcs, spread across 
the image and between the BCGs. Additional blue faint im- 
ages between images 1.5 and 1.6 (see Figure [9]) may also 
be related to this system. We did not find any record of 
past analysis nor redshift information and we assume the 
redshift of this system to he Zs ^ 2 corresponding to an 
equivalent Einstein radius of 28 =t 2" enclosing a mass of 
2.07 =b 0.10 X lO^^M©. This also corresponds to the higher 
redshift assumed for system 2, which is estimated to be a 
drop-out galaxy at Zs 3.5 helping us to constrain the 
fit. See Figure [9] for more details. The reference centre of 
our analysis is fixed between the two cD galaxies at: RA = 
06:47:50.23, Dec = +70:14:55.37 (J2000.0), where one arc- 
second corresponds to 6.67 kpc at the redshift of this cluster. 
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Figure 8. 2D surface mass distribution (k), in units of the crit- 
ical density, of MACS J0454. 1-0300. Contours are shown in hn- 
ear units, derived from the mass model constrained using the 
multiply- lensed images seen in Figure |7| Axes are in ACS pixels 
(0.05" /pixel), and a 20" scale bar is overplotted. 



Figure 10. 2D surface mass distribution (k,), in units of the crit- 
ical density, of MACS J0647. 7+7015. Contours are shown in lin- 
ear units, derived from the mass model constrained using the 
multiply-lensed images seen in Figure |9] Axes are in ACS pixels 
(0.05" /pixel), and a 20^^ scale bar is overplotted. 




Figure 9. Galaxy cluster MACS J0647.7+7015 (z = 0.591) im- 
aged with Hubble/ACS F555W and F814W bands. The critical 
curves overlaid in white correspond to the six images we identi- 
fied as belonging to system 1, at an estimated source redshift of 
Zs ~ 2, enclosing a critical area with an effective Einstein radius 
of ~ 190 kpc, at the redshift of this cluster. The second system 
consists of 3 drop-out candidate images, corresponding to a red- 
shift of Zs ~ 3.5. 



3.6 MACS J0717.5+3745 

The galaxy cluster MACS J0717.5+3745 (z = 0.55) was 
analysed recently by Zitrin et al. (2009a). This very X-ray 
luminous cluster is the denser north-western region of the 
large-scale filament found by Ebeling et al. (2004) and it is 
the largest known lens, with an equivalent Einstein radius 
of 55 ''and a mass of 7.4 ±0.5 x lO^^M©. 13 multiply-lensed 
systems were used to constrain the fit as seen in Figure [TT] 
taken from Zitrin et al. (2009a). Due to its very large critical 
area, this cluster is a great target for finding high-z objects, 



and has been recently proposed for uncovering early stars 
(or "dark stars", Zackrisson et al. 2010). 



3.7 MACS J0744.8+3927 

The luminous X-ray cluster MACS J0744.8+3927 (z = 
0.698) is the highest-^ cluster of this sample, comprising 
several obvious close pairs of multiply-lensed galaxies which 
are immediately visible throughout the frame, with which 
we begin our modelling process (see Figure 13). Systems 2 
and 4 are likely drop-outs and we estimate their redshift 
as Zs ^ 3.5 which helps to constrain our model. These 
systems correspond to similar lensing distance ratios thus 
basing this assumption. We did not find record of past 
analysis yet recently a spectroscopic redshift of a resolved 
galaxy behind this cluster was published (Jones et al. 2009; 
Zs = 2.21, marked as "SI" in Figure 13 here). Unfortu- 



nately, no counter images were reported. As can be seen in 
Figure this arc lies within the critical curve for a source 
at Zs — 2.2 according to our model, and therefore should be 
multiply-lensed and we mark other possible counter images 
of this galaxy. Another option, which would favour a steeper 
model, is that the three clumps seen in this arc are counter 
images of the same galaxy, though we find this option to 
be less probable according to the velocity map presented 
in Jones et al. (2009), which shows that the three clumps 
have different velocities. In addition, another less probable 
case is that the critical curve for a source at Zs — 2.2 is 
much smaller, though this does not agree with the other 
assumed redshifts for the outer blue arcs (system 1 for ex- 
ample), which are not likely be at a higher redshift than 
Zs ^ 3 since they are closer to the cluster center than sys- 
tem 2 but following the same symmetry, yet clearly seen in 
the F555W band. Deeper imaging in more passbands will 
help to uniquely solve the strong-lensing in this cluster. We 
derive an equivalent Einstein radius of 31 di 2'' for a source 
redshift of Zs — 2.2, enclosing a mass of 3.1±0.1 X lO^^M©. 
The reference centre of our analysis is fixed on the cD galaxy 
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Figure 11. Galaxy cluster MACS J0717.5+3745 {z = 0.546) im- 
aged with Hubble/ACS F555W, F606W, and F814W bands. The 
34 multiply lensed images identified by our model are numbered 
here. The white curve overlaid shows the tangential critical curve 
corresponding to the distance of system 1 at an estimated red- 
shift of z ~ 2 — 2.5, and which passes through several close pairs 
of lensed images in this system. The larger critical curve over- 
laid in red corresponds to the larger source distance for the red 
dropout galaxy number 5, at the estimated photometric redshift 
of 2; ~ 3.5 — 4. This large tangential critical curve encloses a very 
large lensed region equivalent to ~ 400 kpc in radius at the red- 
shift of the cluster. Figure was originally published in Zitrin et 
al. 2009a. 
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Figure 12. 2D surface mass distribution (k), in units of the crit- 
ical density, of the inner central region of MACS J0717. 5+3745. 
Contours are shown in linear units, derived from the mass model 
constrained using 34 multiply-lensed images seen in Figure pT] 
Axes are in ACS pixels {0.05" /pixel), and a 20^^ scale bar is over- 
plotted. Note that the central mass distribution of is rather fiat 
refiecting the unrelaxed appearance of this cluster. Figure was 
originally published in Zitrin et al. 2009a. 




Figure 13. Galaxy cluster MACS J0744.8+3927 (z = 0.698) im- 
aged with Hubble/ACS F555W and F814W bands. The critical 
curves overlaid in red correspond to systems 2 and 4, at an esti- 
mated source redshift of Zs ~ 3.5. The inner critical curve overlaid 
in white corresponds to a source redshift of Zs — 2.2 (arc "SI"; 
measured spectroscopically to be at this redshift by Jones et al. 
2009). We mark possible counter images of this arc and other 
multiply-lensed images we uncovered throughout the frame. 
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Figure 14. 2D surface mass distribution (^c), in units of the crit- 
ical density, of MACS J0744. 8+3927. Contours are shown in hn- 
ear units, derived from the mass model constrained using the 
multiply-lensed images seen in Figure [Ts] Axes are in ACS pixels 
(0.05'' /pixel), and a 20^^ scale bar is overplotted. 

at: RA = 07:44:52.80, Dec = +39:27:26.50 (J2000.0). One 
arcsecond corresponds to 7.17 kpc at the redshift of this 
cluster. 

3.8 MACS J0911. 2+1746 

In the galaxy cluster MACS J0911.2+1746 (z = 0.505) not 
many prominent arcs are seen, in agreement with the small 
critical area we derived (Oe = lllf', for a source redshift 
of Zs ^ 2). We found no record of previous SL analysis of 
this cluster nor arcs redshift information. Only a few arcs are 
seen throughout the frame, from which we are able to match 
7 multiply-lensed images which are iteratively incorporated 
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Figure 15. Galaxy cluster MACS J0911. 2+1746 {z = 0.505) 
imaged with Hubble/ACS F555W and F814W bands. The 7 
multiply-lensed images identified by our model are numbered 
here. The white curve overlaid shows the tangential critical curve 
corresponding to the distance of system 1 at an estimated red- 
shift of 2; ~ 2. This cluster comprises the smallest Einstein radius 
and mass of the analysed sample. The area enclosed within the 
marked critical curve corresponds to an equivalent Einstein radius 
of 68 kpc at the cluster redshift. 
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Figure 16. 2D surface mass distribution (/^), in units of the crit- 
ical density, of MACS J0911. 2+1746. Contours are shown in lin- 
ear units, derived from the mass model constrained using the 
multiply-lensed images seen in Figure^] Axes are in ACS pixels 
(0.05" /pixel), and a 20^^ scale bar is overplotted. 



into the model, belonging to 3 distant galaxies. Two of these 
galaxies (systems 1 and 2) are to our estimation at Zs ^ 2. 
We find that the mass enclosed within the critical curves 
(for a source redshift of Zs - 2) is 0.28t° °? x lO^^M©. This 
critical curve and the enclosed mass are both the smallest of 
the current sample. See Figures [15] and [16] for more details. 
The reference centre of our analysis is fixed near the centre 
of the ACS frame at: RA = 09:11:10.30, Dec = +17:46:39.33 
(J2000.0), where one arcsecond corresponds to 6.17 kpc at 
the redshift of this cluster. 




3.9 MACS J1149.5+2223 

This cluster was first analysed by us (see Zitrin & Broad- 
hurst 2009) noting several large blue spiral galaxy images 
which are clearly visible near the central brightest cluster 



galaxy (Figure 17). Shortly after, spectroscopic redshifts 



were published (Smith et al. 2009) in full agreement with 
the predictions of our analysis: we assumed ^ :^ 1.5 for 
the spiral-galaxy (system 1), and z ~ 2 for the outer blue 
images (system 3 in Figure 17), which were later verified 
to be at z = 1.49 and z = 1.89, respectively. Many other 
faint lensed galaxies are also visible, most of which we have 
been able to securely identify as belonging to 10 sets of 
multiply-lensed background galaxies (Figure 17); see Zitrin 



& Broadhurst (2009) for more details. We find that the crit- 
ical curves for a source redshift of z :^ 2 enclose a mass of 
1.71=b0.20x IO^^Mq, and have an equivalent Einstein radius 
of 27 ±3''. 



Figure 17. Large scale view of the multiply-lensed galaxies iden- 
tified by our model in MACS J1149. 5+2223 (z = 0.544). In addi- 
tion to the large spiral galaxy system 1, many other fainter sets of 
multiply lensed galaxies are uncovered by our model. The white 
curve overlaid shows the tangential critical curve corresponding 
to the lensing distance of system 1. The larger critical curve over- 
laid in blue corresponds to the average distance of the fainter 
systems, passing through close pairs of lensed images in systems 
2 and 3. This large scale elongated "Einstein ring" encloses a 
very large critically lensed region equivalent to 170 kpc in radius. 
For this cluster one arcsecond corresponds to 6.4 kpc, with the 
standard cosmology. Figure was originally published in Zitrin &; 
Broadhurst 2009. 
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Figure 18. 2D surface mass distribution (k), in units of the criti- 
cal density, of MACS J1149. 5+2223. Contours are shown in hnear 
units, derived from the mass model constrained using 33 multiply 
lensed images seen in Figure pT| Note, the central mass distribu- 
tion is shallow, and rounder in shape than the critical curves. Fig- 
ure was originally published in Zitrin &; Broadhurst 2009. Axes 
are in ACS pixels (0.05" /pixel), and a 20^^ scale bar is overplot- 
ted. 

3.10 MACS J1423.8+2404 

The galaxy cluster MACS J1423.8+2404 (z = 0.543) was 
recently analysed by Limousin et al. (2009), finding 3 sets 
of multiply- lensed galaxies. System 1 here was spectroscop- 
ically measured by them to be at Zs = 2.84 and system 
2 here was spectroscopically measured by them to be at 
Zs = 1.78. We incorporate these images in order to fully 
constrain our model, with an additional locally-lensed arc 
found here (system 3). The corresponding critical curves 
for a source redshift Zs ^ 2 are overlaid in Figure [l9] 
enclosing a mass of 1.3 =b 0.40 x 10^^ Mq and yielding an 
effective Einstein radius of = 20 =b 2^^ This Einstein ra- 
dius is in agreement with that quoted by Limousin et al. 
(2009), yet the mass is lower than quoted by them, but fully 
agrees with the Mass/Critical- area relation seen in Figure 
[27l which supports our measurement. The reference centre 
of our analysis is fixed near the centre of the ACS frame 
at: RA = 14:23:48.05, Dec = +24:05:00.23 (J2000.0), where 
one arcsecond corresponds to 6.40 kpc at the redshift of this 
cluster. 




Figure 19. Galaxy cluster MACS J1423.8+2404 (z = 0.543) im- 
aged with Hubble/ACS F555W and F814W bands. The multiple- 
images used in our model are marked on the image. System 1 was 
spectroscopically measured by Limousin et al. (2009) to be at 
Zs = 2.84 and system 2 was spectroscopically measured by them 
to be at Zs = 1.78. On the image we overlay the critical curves for 
a source redshift of ~ 2, enclosing a critical area with equiva- 
lent Einstein radius of ~ 130 kpc at the cluster redshift. 
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3.11 MACS J2129.4-0741 

In the galaxy cluster MACS J2129.4-0741 (z = 0.589) var- 
ious arcs are seen throughout the image. One spectacular 
system consists of 6 lensed images in similar colours as the 



cluster members (system 1; see Figure 21). We use this sys- 
tem and the large arc (system 2) further away from the cen- 
tre to fully constrain the model. Estimating that the outer 
blue arc (system 2) is at 2— 2.5, correspondingly the six- 
times lensed galaxy of system 1 is at a redshift of Zs 1 — 1.5. 
The outer critical curve, corresponding to system 2 encloses 
an area of an effective radius of 6e = 37 d= 2'' and a mass 
of 3.4l°;^ X lO^^M©. The reference centre of our analysis is 
fixed at: RA = 21:29:26.123, Dec = -07:41:27.28 (J2000.0), 



Figure 20. 2D surface mass distribution (ac), in units of the crit- 
ical density, of MACS J 1423.8+2404. Contours are shown in hn- 
ear units, derived from the mass model constrained using the 
multiply- lensed images seen in Figure [19] Axes are in ACS pixels 
(0.05'' /pixel), and a 20^^ scale bar is overplotted. 



where one arcsecond corresponds to 6.66 kpc at the redshift 
of this cluster. 



3.12 MACS J2214.9-1359 

The galaxy cluster MACS J2214.9-1359 (z = 0.503) com- 
prises several prominent large blue arcs, 4 of which we match 
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Figure 21. Galaxy cluster MACS J2129. 4-0741 (z = 0.589) im- 
aged with Hubble/ACS F555W and F814W bands. The large 
white critical curve corresponds to system 2, at an estimated red- 
shift of Zs ~ 2, enclosing a critical area of an effective Einstein 
radius of ~ 250 kpc at the redshift of this cluster. Comprising 6 
remarkable images in the centre of the image, system 1 is at a 
redshift of 2:5 ~ 1, whose critical curve is overlaid in blue. 




1000 1500 
Pixels 



Figure 22. 2D surface mass distribution (k), in units of the crit- 
ical density, of MACS J2129. 4-0741. Contours are shown in lin- 
ear units, derived from the mass model constrained using the 
multiply- lensed images seen in Figure [21] Axes are in ACS pixels 
{0.05" /pixel), and a 20^^ scale bar is overplotted. 



as system 1 (see Figures 23 and 35). This system might in- 
clude an additional image, with mirror symmetry to 1.4, if 
the nearby foreground object is prominently included, yet 
we do not use this image for constraining our model. We 
estimate the redshift of systems 1 and 2 as ~ 2, which 
encloses a mass of 1.25±O.lOxlO^^M0 in a critical area with 
an effective Einstein radius of 6e — 23 ^ 2!' . The reference 
centre of our analysis is fixed at: RA = 22:14:56.59, Dec = 
-14:00:17.23 (J2000.0), where one arcsecond corresponds to 
6.16 kpc at the redshift of this cluster. 



Figure 23. Galaxy cluster MACS J2214.9-1359 {z = 0.503) im- 
aged with Hubble/ACS F555W and F814W bands. The critical 
curve overlaid in white corresponds to system 1, enclosing a criti- 
cal area of an effective Einstein radius of ~ 140 kpc at the redshift 
of the cluster, for an estimated source redshift of Zg ~ 2. 
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Figure 24. 2D surface mass distribution (/^), in units of the crit- 
ical density, of MACS J2214. 9-1359. Contours are shown in lin- 
ear units, derived from the mass model constrained using the 
multiply- lensed images seen in Figure [23] Axes are in ACS pixels 
{0.05'' /pixel), and a 20^^ scale bar is overplotted. 



4 DISCUSSION 

Here we examine the SL properties of the whole sample. The 
Einstein radii, enclosed mass, and high-z magnifications are 
listed in Table |2] and discussed in the following subsections. 



4.1 The Einstein Radius Distribution 

The 12 clusters have effective Einstein radii in the range, 
11^^ to 55^^ with a median value of 27.5" , (and a mean value 
of 28 =b 3.6^^) assuming a fixed source redshift of Zs — 2, 
see Figure [25] This corresponds to a range of physical radii, 
68 kpc to 353 kpc, with a median value of 180 kpc, when 
transforming to the redshift of each lens. In the hierarchi- 
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Table 2. Summary and results of the SL analysis. Part of the following data are based on Ebeling et al. (2007; see also TableQ. Column 
3: Effective Einstein radius, in arcseconds. Simply the square root of the area enclosed within the critical curves divided by tt; Column 
4: Effective Einstein radius, in kpc, for Zs ~ 2. An uncertainty of Az zb 0.5 in the estimated source redshift results in a typical ~ 10% 
deviation in the lensing distance ratio, and usually up to such deviation also in the Einstein radii. Column 5: Mass enclosed within the 
critical curves, in lO^^M©. The quoted errors correspond to the different models. They do not include the source redshift uncertainty 
where exists. An uncertainty of Az ~0.5 can cause a typical mass uncertainty of ^^q%; Column 6: M/Lb-, in {M/L)q, fluxes were 
converted to luminosities using the LRG template described in Benftez et al. (2009); Column 7: Lower limit of the highly magnified 
(> xlO) area for a high- 2; source at Zs ~ 8, in square arcminutes; Column 8: image-plane RMS, in arcseconds; for more details see also 
Table [1] 
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cal model large-scale perturbations collapse recently and 
thus should be found relatively locally. This sample includes 
very large Einstein radii, exceeding even the most impressive 
lenses previously studied at lower redshifts as analysed by 
Broadhurst & Barkana (2008). More recently, Richard et al. 
(2009b) carefully studied the lensing properties of 20 mainly 
undisturbed clusters at lower- 2;, selected to have SL features 
in Hubble data and measured X-ray data, for which their 
mean Einstein radius is Oe — 14. 45^^ a factor of ^two smaller 
than our sample. This difference reflects the larger masses of 
the MACS sub-sample that we studied, which is purely X- 
ray flux selected and restricted to z > 0.5. This MACS sam- 
ple therefore does not suffer from a lensing-selection bias, 
but the effect of projection bias must be taken into account 
when evaluating 2D lensing observations with 3D mass pro- 
files predicted by theory (e.g., Hennawi et al. 2007, Oguri & 
Blandford 2009, Sereno, Jetzer & Lubini 2010), as discussed 
in ^4.41 below. 



Most of these clusters have resisted analysis by strong 
lensing despite the long availability of the Hubble data. Our 
success in "cracking" these irregular lenses and defining their 
critical curves is encouraging for the application of our ap- 
proach to complex unrelaxed clusters more generally. For 
such clusters, where the deflection fields are not symmetric, 
the identification of multiple images must be aided by mod- 
els that allow for flexibility in describing the mass distribu- 
tion, as described in section |3] free of the strictures imposed 
by the use of idealised elliptical potentials in conventional 
models. 



Theor. Dist. 
- - - Obs. Median 
Theor. Median 




Figure 25. Einstein Radius distribution. This histogram shows 
the number of clusters per bin, where 6 equally-spaced bins were 
used to divide the range 11^^ to 55^^ The solid curve is the ex- 
pected distribution of Einstein radii as calculate in g4.4| for the 
ACDM model, incorporating both the scatter on the Lx — M 
relation and the projection bias from lensing. We multiply the 
resulting theoretical dN/dOe curve by the width of the bins to 
normalise it. The observed distribution is skewed to larger Ein- 
stein radii than predicted. The median values of both the observed 
and the theoretical distributions are plotted on the histogram in 
a blue dashed line, and a black dash-dotted line, respectively. 



4.2 Central Mass Distributions 

The masses enclosed by the critical curves range from 

median value 

14 7 



2.8: 



X 10''' Mo to 7.4±0.5 X IO'^Mq with 



of 1.9 X lO'^Mo, and a mean value of 2.4l°;^ x lO'^'M© as 
shown in table [2] These are calculated for each cluster by 
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Figure 26. Comparison of the cumulative distribution of ob- 
served Einstein Radii (blue solid stairs) and the theoretical dis- 
tribution predicted by the ACDM model (black solid curve) de- 
scribed in §4.4| The median values of both the observed and the 
theoretical cumulative distributions are plotted on the histogram 
in a blue dashed line, and a black dash-dotted line, respectively. 
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Figure 27. Mass enclosed within the critical curves as a function 
of the effective Einstein radii. The theoretical relation for sym- 
metric mass distribution, M oc r^, is tightly followed by the data, 
indicating that in general the central mass distributions of these 
clusters are well behaved. 



integrating the surface mass density distribution within the 
2D model critical curve, scaled to a source redshift of Zs ^ 2. 
We examine the relation between these Einstein masses, 
Mein enclosed within the critical curves, and the effective 
Einstein radii derived above. Theoretically this should of 
course simply scale as Mein oc Ol for symmetric lenses. Here 
the lenses are not symmetric, but as can be seen in Figure 
|27|the quadratic relation is tightly followed, indicating that 
asymmetry is not very significant in terms of the mass distri- 
bution, and furthermore the measured or assumed redshifts 
of z ^ 2 for the relevant systems is a reliable estimate. In 
fact, we see clearly in that the 2D mass distributions are 
in general noticeably rounder than the critical curves which 
are very sensitive to substructure, as can be seen by com- 
paring the 2D mass distributions with the critical curves. 




Figure 28. Mass-to-light ratio enclosed within the critical curves 
as a function of the Einstein radius (in arcseconds). A general 
trend of increasing M/L with the system scale size is apparent. 




Figure 29. Reproduction of system 1 in MACS J0018. 5+1626 
by our model. Since the arc is too faint to be noticed here we lens 
the middle arclet and paint in red the triplet outcome. 



shown in Figures [T] - 124| Further encouragement for the ac- 
curacy of SL mass estimates is presented in recent work by 
Meneghetti et al. (2010), showing that SL methods based on 
parametric modelling are accurate at the level of few percent 
at predicting the projected inner mass. 

Richard et al. (2009b) found a mean mass of 1.95 x 
lO^^M© enclosed mass (within R<250 kpc) for their sam- 
ple of 20 clusters, while Smith et al. (2005) examined the 
mass distribution of 10 X-ray selected clusters, and found a 
similar mean mass of 1.86 x lO^^M© (see also Richard et al. 
2009b). These values are in good agreement with our sam- 
ple median mass (1.9 x lO^'^M©), and slightly lower than 
our sample mean mass which is boosted by few extremely 
massive clusters (e.g., MACS J0717.5+3745). 

We examine the behaviour of the central mass-to-light 
ratio, M/Lb, with the Einstein radius and shown in Fig- 
ure [28] The luminosity is summed over all cluster sequence 
galaxies identified as described in section |3] A clear corre- 
lation is obtained such that M/Lb increases with the Ein- 
stein radius, and with values quite typical of other massive 
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Figure 30. Reproduction of system 1 in MACS J0025. 4-1222 by 
our model, by delensing image 1.1 into the source plane, and re- 
lensing the source-plane pixels onto the image plane to accurately 
form the other images. 

clusters (Sarazin 1988, Rines et al. 2004, Medezinski et al. 
2010). This scaling towards higher M/Lb provides confi- 
dence in our method. At the high mass end our result is 
similar to the peak M/L found in the detailed weak lensing 
profile studies of Medezinski et al. (2010). Measuring M/Lb 
in high-z clusters is of interest also for the added insight on 
the amount of DM initially associated with individual galax- 
ies as compared with the overall cluster DM component (see 
also Sarazin 1988). 

4.3 The Magnification Distribution 

As can be seen in Table [2] most of the sample clusters pro- 
vide large regions (^ 2.5n^) of highly magnified (> xlO) 
sky, useful for the search of the first stars and galaxies (e.g., 
Zackrisson et al. 2010). One has to be cautious when mak- 
ing such a statement since unlike the critical curves and the 
enclosed mass, the magnification is sensitive to the mass 
profile which requires source redshift estimation for several 
sources, which we have here for only a few clusters (section 
|3|. Still, as explained in §2, reasonable constraints are put 
on the mass profile slope by the image-plane minimisation, 
which generally has a broad minimum at the same point in 
the parameter space where the slope is also approximately 
correct. In addition, we are able to put a lower boundary 
on the high magnification area, by choosing the steeper pro- 
files which naturally generate lower magnifications, or by 
assuming a boosted redshift for the main multiply-lensed 
system decreasing the lensing distance ratio for high-z galax- 
ies. Thus we constrain the area of high magnification for 
a high-z source at the current limit of ^ 8, showing 
that most of this sample clusters are excellent targets for 
high- 2; galaxy searches due to these large high-magnification 
areas and lens power. We estimate, based on previously 
analysed clusters (e.g., CI 0024+1654, Zitrin et al. 2009b; 
MACS J1149.5+2223, Zitrin & Broadhurst 2009; MACS 
J0717.5+3745, Zitrin et al. 2009a) that in practice the ar- 
eas of high magnification are about 10-20% higher than the 
lower limits presented in Table |2] with an upper limit of 
about 40-50%. 

4.4 Comparison with ACDM 

We generally follow Broadhurst & Barkana (2008; see also 
Sadeh & Rephaeli 2008) in constructing the theoretically 




Figure 31. Reproduction of system 1 in MACS J0257. 1-2325 
by our model, by delensing image 1.1 into the source plane, and 
relensing the source-plane pixels onto the image plane to accu- 
rately form the other images. Note, our model predicts an extra 
tiny image near the core of the cD covered by its light. 
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Figure 32. Reproduction of system 1 in MACS J0647.7+7015 by 
our model, by delensing image 1.1 into the source plane, and re- 
lensing the source-plane pixels onto the image plane to accurately 
form the other images. 

predicted distribution of Einstein radii. In particular, we 
adopt the NFW parameters measured by Neto et al. (2007) 
for simulated halos, although we reduce cat by 10% accord- 
ing to the results of Duffy et al. (2008) who used the more 
recently measured values of the cosmological parameters. 
Studies based on large numerical simulations (Zhao et al. 
2003,2009, Gao et al. 2008) have found for massive halos a 
rather weak decline of cat with increasing redshift; we neglect 
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Figure 33. Reproduction of system 1 in MACS J0717.5+3745 
by our model. Published originally in Zitrin et al. 2009a. 
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Figure 34. Reproduction of system 1 in MACS J1149. 5+2223 
by our model. Published originally in Zitrin &: Broadhurst 2009. 

this decline, which renders our results (in terms of the dis- 
agreement with observations) conservative. We then correct 
the concentration parameter distribution using Hennawi et 
al. (2007) to obtain the effective parameters for the popula- 
tion of clusters observed in projection. Note though that here 
we are considering a sample selected by X-ray flux (not lens- 
ing cross-section), so there is no le using bias, but there still 
is the projection bias; this refers to the fact that when halos 
are seen in projection (at randomly distributed angles), a 
different effective projected cat is measured than would be 
obtained from a spherical analysis of the halo density pro- 
file. In particular, projection increases the mean cn and adds 
substantially to the scatter compared to the distribution of 
the spherically-measured cn- Neto et al. (2007) split their 
halo sample into two groups ("relaxed" and "unrelaxed") 
and analysed the statistics of each group separately, so we 
do the same, and in the end combine the two groups ac- 
cording to their relative numbers in the simulation by Neto 
et al. (2007), obtaining our predicted results for the total, 
combined population of halos. This can then be compared 
with the observed sample, which does not select for relaxed 
clusters but includes both relaxed and unrelaxed ones. 

In order to approximately simulate a flux-selected sam- 
ple, we need to convert Lx to halo mass. We use the power- 
law relation from Reiprich & Bohringer (1999; with an up- 
dated Hubble constant), along with the observed scatter of 
Lx for a given mass, which we model as a lognormal distri- 
bution with a typical scatter that corresponds to a factor of 
1.5 . Note that this observed relation is based on M500, while 
the results from the above simulations refer to M200, so we 
must include the relation between M500 and M200 which is 
itself a function of ca^. An important caveat is that the ob- 
served Lx — M relation was measured for relaxed clusters 
(whose mass could be estimated from the X-ray emission as- 
suming hydrostatic equilibrium). We expect that unrelaxed 
clusters would tend to have an unusually high Lx for a given 
mass (e.g., in a post-merger phase), which would insert more 
low-mass clusters within the sample than we assume, leading 
to smaller predicted Einstein angles, increasing the discrep- 
ancy between the theory and the observations (again making 
our results conservative). 

On the other hand, it should be noted that much of 
the scatter in the Lx — M relation probably arises from 
cool cores in the center (since it is signiflcantly reduced if 



the inner core is excised), implying a possible bias towards 
relaxed and more concentrated systems. However, such a 
bias seems less probable in our case due to the high X-ray 
luminosity, the relatively high redshift, and the clear spread- 
out appearance of most of the clusters analysed here. 

The observed properties of each halo (flux, Einstein ra- 
dius distribution, etc.) change slowly with redshift, so to 
simplify the calculation we approximate the clusters as all 
lying at a typical redshift z = 0.55. However, the halo abun- 
dance changes rapidly with redshift, so we calculate the in- 
tegrated halo mass function for all halos observed at 2; > 0.5, 
in the fraction of the sky corresponding to the MACS survey 
(around a quarter of the sky). We use the Sheth & Tormen 
(1999) formula, which accurately flts the halo mass function 
(i.e., the distribution of M200) measured in simulations. 

This halo mass function was convolved with with the cn 
distribution to obtain the distribution of M500, and then add 
the Lx — M500 scatter in order to obtain the halo mass func- 
tion corresponding to selection by a given minimum flux. 
Since the Lx — M500 relation is not strictly valid for unre- 
laxed clusters, we only use the power- law and the scatter 
from this relation but allow some flexibility in the normali- 
sation, choosing an effective minimum flux that corresponds 
to a predicted abundance of 12 clusters (as in the observed 
sample). This is reasonable also since we are interested here 
in testing the predicted Oe distribution, not in testing the 
predicted cluster abundance. Note that the flux selection 
does not affect the distribution of the very highest Oe val- 
ues, since all clusters within the survey solid angle at 2; > 0.5 
that can produce the largest ^e's will be well above the flux 
threshold. 

Given our predicted. X-ray selected mass function in 
terms of M500, we must convolve it with the Oe distribution 
of each halo of mass M500 • From the simulations we have the 
Cn distribution (and thus also the Oe and M500 distribution) 
for a given M200, but we can invert this conditional prob- 
ability using Bayes' theorem. Note that the M200 to M500 
conversion should not really include the projection bias that 
is included in ^e, but having two different cn distributions 
for each M200 would makes things far more complicated. Our 
simplification of using a single cn distribution is reasonable 
since the projection bias typically affects M500 (for a given 
M200) only by 15%. Also, since cn affects both M500 and Oe 
in the same direction, the effect on our flnal result (which 
derives Oe for a given M500) is much smaller. 

Our predicted Oe distribution is compared to the ob- 
servations in Figure [26] We compare the two cumulative 
distributions A^(> ^e), i-e., the total number of clusters ex- 
pected above Oe , where the theoretical distribution has been 
normalised (as noted above) to agree with the total of 12 
clusters in the sample. The observed distribution clearly lies 
at higher Einstein radii compared with the predicted distri- 
bution. One indication of this is the most discrepant cluster, 
J071 7. 5+3745 with Oe = 55''. The theoretical calculation 
yields a probability of only 3.4% of flnding such a large Oe 
in this z > 0.5 cluster sample, giving a signiflcant (though 
not extreme) 2-cr discrepancy. More generally, the observed 
distribution resembles the predicted one except offset to- 
wards higher angles by about a factor of 1.4 (comparing, 
for example, the two medians). The standard K-S statistic 
for comparing two distributions gives a probability of 7.4% 
that the observed distribution is drawn from the predicted 
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one. This is also around a 2-a discrepancy, largely indepen- 
dent of the above number since the K-S statistic focuses on 
the middle portion of the probability distributions and is 
insensitive to the edges. The relatively low significance of 
the K-S discrepancy is an inevitable result of the Poisson 
fluctuations expected with a sample of only 12 clusters. 

The discrepancy that we find, of Einstein radii that are 
too high by about a factor of 1.4, appears smaller than 
the roughly factor of 2 difference found by Broadhurst & 
Barkana (2008). This may be explained at least partially by 
the difference between a lensing selected and an X-ray se- 
lected sample. A more robust comparison would be possible 
with an X-ray selected cluster sample in which the clusters 
have reliable virial masses measured using weak lensing. An- 
other caveat is that we have compared the observed Einstein 
radii with the predictions based on pure dark matter sim- 
ulations. The possible effect of baryons on the halo density 
profile has been difficult to study with simulations, since 
hydrodynamic simulations that produce a large effect do 
so along with a central baryon concentration that disagrees 
with observations and is due to the "overcooling" problem 
in simulations. Analytical models suggest that even without 
producing a central baryon concentration, the complex cou- 
pled history of baryon and dark matter accretion may allow 
the baryons to significantly affect the final central density 
profile, potentially alleviating the discrepancy at least par- 
tially (Barkana & Loeb 2010). On the other hand, recent 
simulations suggest that the discrepancy actually increases 
when AGN feedback is taken into account to overcome the 
overcooling problem (Duffy et al. 2010). 



4.5 Effects of Uncertainties on the Results 

In this work we determine the Einstein radii and projected 
masses for the 12 clusters of the z > 0.5 MACS sample, 
and compare these to simulations. We make use of our 
well-established modelling method to identify many sets of 
multiply-lensed images which are in turn used to constrain 
the models, and determine the Einstein radii and projected 
masses for Zs ^ 2. Typical uncertainties in the source red- 
shift estimates of Az ±0.5 impose only minor uncertainties 
on the comparison with ACDM. Firstly, such an uncertainty 
does not enable a unique determination of the profile slope, 
which we do not attempt to fully constrain here since both 
the Einstein radius and the projected mass are indifferent 
to the mass profile slope and are determined by the data. 
Secondly, overestimating a source redshift would in practice 
increase the projected mass and the observed Einstein radius 
for a source at 2;^ = 2, while maintaining the relation seen in 
Figure |27| and thus resulting only in growth of the discrep- 
ancy from ACDM simulations which are not dependent on 
the observed mass. On the other hand, underestimating a 
source redshift would in practice decrease the observed Ein- 
stein radius and projected mass for a source at Zs = 2 by 
less than 10%, thus insignificantly influencing the results. 



5 SUMMARY 

The MACS ^ > 0.5 sample has proven to be of great value 
for many very different studies spanning the X-ray through 



radio spectrum. We have extracted additional useful in- 
formation from this sample by solving the strong lensing 
for these clusters. Previous discoveries from this sample in- 
clude the largest lensed images of a highly- magnified distant 
spiral galaxy (MACS Jl 149.5+2223; Zitrin & Broadhurst 
2009, Smith et al. 2009), the largest known lens (MACS 
J0717. 5+3745; Zitrin et al. 2009a), lensed sub-mm sources 
(Takata et al. 2003, Borys et al. 2004, Berciano Alba et al. 
2007, 2009), and another "bullet cluster" (MACS J0025.4- 
1222; Bradac et al. 2008b). 

We presented mass models and the corresponding crit- 
ical curves for the 12 high- 2; X-ray luminous clusters of the 
sample via strong- lensing analysis in HST/ACS images. Sev- 
eral of these clusters have only a few strongly- lensed im- 
ages and corresponding low masses and Einstein radii, while 
most are very massive and rich with multiply-lensed arcs 
due to large critical area and high lens power. The equiv- 
alent Einstein radii of this sample range from 11^^ to 55^^ 
with a median value of ^ 28" (similar to the mean value), 
or 180 kpc. The corresponding masses enclosed within these 
curves range from 2.8t°;? x lO^^M© to 7.4 ± 0.5 x lO^^M© 
with a median value of 1.9 x lO^^M©, and a mean value of 
2.4+°;^ X lO^^M©. We find that the enclosed mass follows 
tightly the quadratic relation with the equivalent Einstein 
radius, indicating that the deviation from symmetry is not 
prominent in these clusters and that our measurements are 
highly accurate given the measured or assumed redshifts. In 
addition, we have shown that M/Lb increases as expected 
with the scale size of the system, or Einstein radius. 

We compare these results to the predictions of ACDM 
taking into account projection biases and find that these 
predictions fall short of our measured Einstein radii by a 
factor of 1.4. In addition, the standard K-S statistic gives 
a probability of only 7.4% that the observed distribution 
is drawn from the predicted one, corresponding to a ^ 2cr 
discrepancy. It is apparent that the observed Einstein radii 
and implied high masses are not likely at high redshifts in 
the context of the ACDM model when no lensing-selection 
biases are involved, showing again possible tension with the 
prediction of the standard model. Finally, as a result of their 
unrelaxed mass distributions, most of these clusters cover a 
very large area (> 2.5n') of high magnification (> xlO) 
making them primary targets for high-z galaxy searches, for 
which substantial HST time will shortly be forth-coming. 

It should be acknowledged that theoretical predictions 
may clearly entail further uncertainties than those taken into 
account here, and a revised and more thorough analysis is 
needed in order to obtain higher significance results once 
source redshifts are available. Such comparisons should also 
be independently applied to other large samples, so that a 
significant number of clusters are compared to the ACDM 
model in order to statistically support any claimed discrep- 
ancy. 
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Figure 35. Reproduction of system 1 in MACS J2214.9-1359 by 
our model, by delensing image 1.1 into the source plane, and re- 
lensing the source-plane pixels onto the image plane to accurately 
form the other images. 
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